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We propose new searches that exploit the unique signatures of decaying sterile neutrinos with 
masses below Mw at the LHC, where they can be produced in rare decays of Standard Model 
gauge bosons. We show that for few-GeV-scale sterile neutrinos, the LHC experiments can probe 
mixing angles at the level of 10“^ — through powerful searches that look for a prompt lepton in 
association with a displaced lepton jet. For higher-mass sterile neutrinos, i.e. Mn > 15 GeV, Run 
H can explore similarly small mixing angles in prompt multilepton final states. This represents an 
improvement of up to two orders of magnitude in sensitivity to the sterile neutrino production rate. 


I. INTRODUCTION 

Run I of the Large Hadron Collider (LHC) proved a 
success by any measure, leading to the discovery of the 
Standard Model (SM) Higgs boson [Tl[2] after decades of 
pursuit by previous experiments. Moreover, the successes 
of the SM have now been extended by the LHC experi¬ 
ments to energy scales well above a TeV, as evidenced by 
the exquisite agreement between data and theory. How¬ 
ever, there are key outstanding missing pieces in the SM, 
chief among them the origin of neutrino masses [3], the 
identify of dark matter, and the dynamics responsible for 
the baryon asymmetry [4] . Remarkably, a minimal exten¬ 
sion of the SM with three gauge-singlet, “sterile” right- 
handed (RH) neutrinos can resolve all of these problems 
if the three right-handed neutrinos all lie below the weak 
scale. Known as the neutrino minimal SM (z^MSM) [SH^, 
the new sterile states lie within kinematic reach of various 
laboratory probes [8], leading to the exciting possibility 
that existing and upcoming experiments could shed light 
on the particles responsible for neutrino masses [Ms], 
dark matter [sunHni, and baryogenesis [HHIIBHll] 1. 

After electroweak symmetry breaking, the z^MSM pa¬ 
rameters consist of the masses of the RH neutrinos, Nj, 
and a mixing angle between each flavor I of sterile neu¬ 
trino and each flavor a of active neutrino, VaNi- The 
interactions between Nj and SM fields are then com¬ 
pletely determined by the mixing between neutrinos; Nj 
can be produced through any interaction involving SM 
neutrinos, with two examples shown in Fig. The ma¬ 
trix element is the same as for the corresponding process 
for SM neutrinos, with an extra multiplicative factor of 
VaNi- The same interactions also allow N to decay to 
SM states. As we elaborate in Section [TTj we focus on 
the mixing angle with the muon neutrino, V^Ni- 

The current constraints for a single RH neutrino mix¬ 
ing with are shown in Fig. § [S usHsg, with the 
least explored parameter space being Mn > mij. The 
proposed SHiP [39] and DUNE/LBNF experiments [40| 
are expected to greatly improve sensitivity for N below 


^ A recent review of the collider phenomenology of RH neutrinos 
can be found in Ref. [25]. 



FIG. 1: Feynman diagrams showing possible production 
mechanisms of a right-handed neutrino, N, via mixing with 
SM neutrinos in (left) the decay of a SM gauge boson (VF^ ^ 
i^N); (right) the decay of a H meson. 

the charm mass (and, to a lesser extent, the b mass), as 
shown in Fig.[^ Above the b mass, the dominant produc¬ 
tion mechanism for N is via the decay of gauge bosons, 
i^N and Z uN. Experiments such as SHiP 
lack the center-of-mass energy needed to produce on-shell 
gauge bosons, thus experiments at high-energy colliders 
such as LEP [29l |3T] and the LHC [§5l |38] are needed 
to discover the RH neutrinos. LEP was sensitive to N 
production via the process Z ^ uN, but its sensitivity 
was limited to \VaNi\‘^ < 10“^ for Mn < M^. 

The LHC features large luminosity and cross sections 
for N production via the decay of SM vector bosons, 
but the soft decay products oiWjZ could be mimicked 
by various reducible and irreducible SM processes. Eur- 
thermore, for Mn ^ Afvp, the decay products of N are 
collimated and may fail isolation requirements in conven¬ 
tional searches. Therefore, carefully targeted analyses 
are needed to extract the sterile neutrino signatures. We 
perform a comprehensive study of signatures of sterile 
neutrinos from decay at the LHC in the kinematic 
regime Mn < Mw We propose new searches for two 
generic parts of parameter space, namely the case where 
N decays at a displaced vertex, and the case where N de¬ 
cays promptly. In the low-mass regime, Mn ^15 GeV, 
the sterile neutrino is both long-lived and boosted; its 
leptonic decays will typically fail standard lepton isola¬ 
tion, but instead give a distinctive signature — that of a 
displaced lepton jet (see [4TH48] for a discussion of lep¬ 
ton jets in other contexts). In the higher-mass regime, 15 
GeV < Mn ^ Mw^ the sterile neutrino decays promptly 
and to three separately resolved objects. We find here 
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FIG. 2: Constraints (shaded regions bounded by solid lines) 
on the right-handed neutrino mass, Mn, and squared mixing 
angle, [25]. We show the estimated 2cr reach of our 

two proposed searches with ^/s = 13 TeV and 300 fb~ ^: the 
displaced lepton jet search discussed in Section III (blue, 
dott ed), and the prompt trilepton search discussed in Sec¬ 
tion IV (brown, dot-dashed). The simplified model is de¬ 
fined in Section [I^ For comparison, the reach of the proposed 
DUNE/LBNF (purple, long dashed) and SHIP (red, short 
dashed) experiments are also shown. 


that targeted trilepton searches can uncover sterile Ma- 
jorana neutrinos in spite of significant SM backgrounds 
from lepton fakes. The trilepton analysis performance 
is more powerful than the only current LHC search for 
Mat < Mw, a fi^{N ^ analysis from CMS 

[38] . The estimated reach of our proposed searches at 
LHC 13 with 300 fb~^ of integrated luminosity is shown 
in Fig. 

This work differs from existing proposals for sterile 
neutrino searches at the LHC in displaced vertices [49| 
by exploiting the distinct, boosted kinematics of the final 
state for Mat ^ Mvp. In particular, for low-mass dis¬ 
placed N, we advocate looking for the distinct signature 
of a prompt lepton in association with a displaced lep¬ 
ton jet. For higher-mass N, when its decays are prompt, 
our proposals for prompt trilepton searches could cover 
new ground by improving sensitivity to \V^n\‘^ in the 
mi) < Mn < Mw range by up to two orders of mag¬ 
nitude over existing bounds. This complements existing 
proposals for LHC searches in and trilepton fi¬ 

nal states that typically target Mat > Mz [5QH^ (al¬ 
though the reach of searches at the Tevatron 

and LHC for Mat < Mw are examined in [55l |57l [60]). 
In this high-mass scenario, N are produced through off- 
shell gauge bosons and thus their decay products exhibit 
very different kinematics from Mat < Mw- Unlike many 
trilepton searches that are targeted at discovering Dirac 
neutrinos or distinguishing between various RH neutrino 
models [58l [59l [CTl [63H67] . we exploit the lepton-number- 
violating character of the Majorana neutrinos in fully 
leptonic decays to improve sensitivity to N. Parts of the 
parameter space that cannot be discovered at the LHC 


with our proposals can potentially be probed in future 
lepton colliders, such as the FCC-ee [68l[6^ . 

The rest of this article is structured as follows. We first 
describe in Section |n| the simplified model under consid¬ 
eration, along with the primary RH neutrino production 
and decay modes. In Section [Hi] we discuss the displaced 
lepton jet signature of RH neutrinos, and in Section HYi 
the prospects for prompt multilepton searches. Finally, 
we conclude with the discussion in Section [V] 


II. SIMPLIFIED MODEL 


The simplified model we consider is the SM with the 
addition of a single new RH neutrino N state: 

FalaHN +]^MN^N, ( 1 ) 

where Lq, is the SM lepton doublet of flavor a, Fq, is a 
Yukawa coupling, and H is the Higgs field. In more real¬ 
istic models of neutrino masses, such as the z/MSM, there 
are three flavors of Y, and the above Lagrangian is sup¬ 
plemented with appropriate indices for the RH neutrinos. 
However, typically only one flavor of N is produced at a 
time in any given interaction at an experiment, and so 
we can consider them independently in evaluating discov¬ 
ery potential. Similarly, we assume for simplicity that N 
couples only to a single lepton doublet a. 

After electroweak symmetry breaking, the RH neutrino 
acquires a Dirac mass in addition to the Majorana mass 
and mixes with the SM neutrinos. In the limit M ^ 
Fa{H), the mass eigenstates are 


^OL 

Mat 



( 2 ) 

( 3 ) 


This is the conventional see-saw mechanism for neutrino 
masses. 

In the same limit, the mixing angles between the active 
and sterile species have the scaling 

Thus, any coupling Qj, involving the SM neutrino gauge 
eigenstate will, after electroweak symmetry breaking, 
also lead to a coupling to N equal to VaN- For two 
examples, see Fig. 

The see-saw relation, Eq. (§, gives a naive prediction 
for VaN as a function of Mat for fixed M^. However, V^at 
depends linearly on Fq,, while M^^ depends on Fj, and so 
a cancellation between real and complex parts of Fq, can 
give rise to a much larger angle than the naive scaling 
suggests I70]2. Additionally, extended theories such as 


^ In Ref. | 71| . it was found that this cancelation may not be radia- 
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the inverse see-saw mechanism m can give rise to VaN 
much larger than the naiVe see-saw prediction; therefore, 
we treat VaN as a free parameter in the simplified model. 

In this paper, we consider the specific simplified model 
where N mixes predominantly with the muon doublet, 
L^. This is the simplified model constrained in cur¬ 
rent LHC searches for RH neutrinos isa EH], and is 
most amenable for detection at colliders due to the clean, 
muon-rich signatures. While we also advocate for direct 
probes of Vqn and VrN at the LHC, these suffer from 
larger mis-tag rates than muons^, and due to the dif¬ 
ficulties of estimating backgrounds, the optimization of 
these studies should be undertaken by the LHC experi¬ 
mental collaborations. We show the current constraints 
on the simplified model parameter space (Mn^V^n) ia 
Fig. II 

The weakest limits are for Mn ^ rric, and we therefore 
focus on this part of parameter space. We further 
narrow our investigation to the promising Mn < Mw 
range, due to the large production rate of N from gauge 
boson and/or meson decays. 

Production: In our simplified model, N production and 
decay proceed entirely through the mixing with the 
The most easily detected production modes for N are: 

1. In W boson decays, 

2. When kinematically allowed, in heavy-flavor 
hadron decays such as 

In both of these production mechanisms, N is produced 
in association with (quasi-)prompt muons and other ob¬ 
jects that can be exploited to discriminate signal events 
from SM backgrounds. We show these production modes 
in Fig. [2 In our analyses, we focus on N production 
from W” decay, but comment on the prospects for pro¬ 
duction from decay in ATLAS and CMS in Section[v| 

Decay: N decays via off-shell W^/Z bosons to other SM 
states. For the masses we consider, rric < Mn < Mw^ the 
RH neutrinos have masses sufficiently above the Quan¬ 
tum Chromodynamics (QCD) confinement scale that 
their decay rates can be computed in terms of quark and 
lepton final states. The proper lifetime, ctat, as a func¬ 
tion of Mn and V/at is shown in Fig. Both prompt 


tively stable in the minimal theory for ^ GeV, although the¬ 
ories with additional lepton number symmetries m may avoid 
this problem. 

^ The sensitivity to VgV of current and upcoming neutrinoless 
double-beta decay experiments ED is also expected to be com¬ 
petitive with, if not better than, the reach of high-energy colliders 
such as the LHC. 

We do not consider the sub-dominant process Z —)■ uN due to 
the lower rate and the lack of a prompt lepton for passing the 
trigger. However, this process can become relevant at a future 

e+e“ collider [Ml ESI- 



Mn (GeV) 

FIG. 3: Contours of fixed right-handed muon-neutrino proper 
decay distance, ctn, as a function of mixing angle and mass. 
The contours are: (right, dashed) ctn = 0.1 mm; (center, 
dotted) ctn = 10 cm; (left, dot-dashed) ctn = 10 m. The 
shaded regions are excluded. 


and displaced decays are expected for RH neutrinos, de¬ 
pending on the mass and mixing angle. 

The most promising decays for the detection of N are: 

1. Leptonic decays, N and N 

In this decay mode, N gives a distinctive final state 
of multiple charged leptons, potentially of different 
flavors, and missing energy; 

2. Semileptonic decays, N /i^qq'. In this case, N 

decays are fully visible, allowing a reconstruction 
of Mn from the decay products, but the hadronic 
backgrounds are potentially larger. For Mn < 
Mw^ there is also typically insufficient transverse 
momentum in each object to separately reconstruct 
each of the quarks as jets as well as the charged lep¬ 
ton at the LHC. 


We find that searches for the clean, fully leptonic signa¬ 
tures to be very sensitive to N production at the LHC for 
Mn < Mw‘ We consider only decays to electrons and 
muons. Neglecting charged lepton masses and assuming 
Mn ^ Afw, the partial widths for leptonic decays rele¬ 
vant for our analyses are [8] 


r(iv ^ 


1927r3 

7687r3 


(a 7^ 0 ), 


X (l + 4 sin^ 6*w + 8 sin^ 0w) 


( 5 ) 

(6) 


where Gp is the Fermi constant and O^v is the Weinberg 
angle. While we show this result for illustrative purposes, 
in Sections in and m we use the full mass-dependent 
widths computed numerically in Madgraph 5 m- 


Search Strategies: In this article, we propose two con¬ 
crete search strategies to discover N, depending on the 
mass and mixing angle of the RH sterile neutrino. 
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1. Displaced lepton jet: For masses rUc < ^ 

the RH neutrino is boosted in the decay 
/a^N. Therefore, the decay products of N are col¬ 
limated. In particular, with N the N 

decay products form a lepton jet. The collider sig¬ 
nature is a prompt lepton + a single displaced lep¬ 
ton jet, which we discuss in Section m 


2. Prompt trileptons: For masses Mn ^15 GeV, the 
RH neutrino is typically neither boosted nor sub¬ 
stantially displaced. The best prospects are in the 

final state, which bene¬ 
fits from smaller SM backgrounds due to the lack 
of an opposite-sign, same-flavor lepton pair. While 
others have considered trilepton probes of very high 
mass RH neutrinos, ^100 GeV [58l [59l ICTl6^ 
[67] . our proposal in Section IV is, to our knowledge, 
the first for the key Mjy < Mw region of parameter 
space, showing trilepton searches to be potentially 
more powerful than the current GMS semileptonic 
searches (381. 


In Section |V| we also briefly consider the prospects for 
discovering N produced in B meson decays for masses 
Mn ^ ^6- 


III. LEPTON JETS FROM RH NEUTRINOS 

This section focuses on the low-mass, ^ Mw 

regime. In this mass range the N will typically decay 
at a length cr > 1 mm for masses above a few GeV but 
less than < 20 GeV, and mixing angles near the current 
limit ^ 10“^, as discussed in Sec. O and shown in 

Fig-i Gonsequently the N decay produ^s will not reg¬ 
ister as prompt final states in the LHG detectors. As long 
as the decay occurs within the fiducial volume of either 
ATLAS or GMS, however, N could still be discovered by 
looking for displaced decay products. 

There are two possibilities to consider, and each moti¬ 
vates a different search strategy, but we focus on the first 
of these (see below). We conservatively consider only the 
case where all of the leptons in the final state are muons, 
as this scenario is less likely to be contaminated by re¬ 
ducible fake SM backgrounds than when an electron is 
present in the final state. 

1. ^ which gives a prompt 

/i & displaced V ^ 2/i + E/t signature; 

2. which gives a prompt 
ja & displaced N /i-f hadronic tracks signature. 

In either scenario above, the prompt fi originates from the 
direct decay of the parent boson, and can be used 
to trigger on the event. Of these two final states, the dis¬ 
placed muon and hadronic tracks signature benefits from 
a higher branching ratio of Br(IF ^ hadrons)/Br(IF ^ 
~ b. However, such a signal region could suf¬ 
fer from potential difficulties, including the possibility 
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FIG. 4: The AR separation between the muons originating 
from the N decays for Mn = 5 GeV (top) and 

Mn = 15 GeV (bottom). The N originates from the decay 
IF^ 


that rare displaced hadron decays, and pile-up-originated 
tracks can mimic the signal. This background is found 
to be negligible with sufficiently high track multiplicity 
and stringent cuts on track px and vertex mass [76] , but 
for Mn ^ Mw^ the signal tracks typically do not have 
enough momentum to pass these cuts. This channel may 
be viable with relaxed kinematic criteria, but we lack 
the tools to properly simulate the backgrounds, as ulti¬ 
mately such a signature would be performed using a data- 
driven estimation of SM backgrounds. Therefore, for the 
low-mass Mn we focus on the first N decay possibility, 
namely the prompt /i+ displaced 2/i + ^x- The prospects 
for this final state were also addressed in Ref. [49| , but we 
propose a different search strategy based on the distinct 
signal kinematics discussed below. 

In the range of N masses and mixing angles giving 
rise to displaced decays, Mn Mw implies that N are 
produced with a significant boost from IF decay. As a 
result, the leptons produced in N decays are typically 
collimated. Fig. shows the AR = Arf + AcjP be¬ 
tween the two muons from the N decay for Mn = 5 GeV 
and Mn = 15 GeV. For Mn = 5 GeV the two final-state 
muons typically appear within less than AR = 0.5 of one 
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FIG. 5: An illustration of the final state strategy advocated 
for the low-mass N regime. The N are produced in the de¬ 
cays of N^ with the subsequent displaced decay 

N The muons from the boosted, displaced N 

will appear as a muon jet /xJ instead of as isolated leptons. 


another. A substantial fraction of leptons from N decay 
for this mass would therefore fail the standard lepton iso¬ 
lation criterion of AR < 0.3 — 0.4. Instead, in analogy 
with the collimated hadrons in a QCD jet, this signature 
suggests a lepton jet (LJ) - namely a reconstructed object 
with a more than one leptonic track concentrated within 
a cone of radius Rq. In this study we only address the 
case of muon-rich final states; thus the object of interest 
is a muon jet 

We use the definition of a muon jet, juJ, from Ref. m 

• Definition: A muon jet is a reconstructed object 
with at least two muons registered in the muon 
spectrometer (MS) within AR < Rq = 0.5 of each 
other. 

Thus the final state that we study in this section is 
that of a prompt muon from the W decay + displaced 
jiJ from N decay. Fig. gives a pictorial representation 
of this signature. While existing lepton jet analyses do 
exist, their signal regions fail to constrain this part of 
parameter space as they require two LJs in the event 
EH, whereas our signal gives one per event. 

Our jiiJ selections are inspired by Ref. The kine¬ 
matic requirements we use to define the signal region for 
the low Mtv displaced lepton jet analysis are related: 


^ Mixed electron-muon jets are also possible, although electrons 
are more easily faked by hadron tracks, so we consider the cleaner 
muon signature to avoid fake backgrounds. 


1. Trigger on a single muon, requiring a prompt muon 
track with px > 24 GeV; 

2. Any remaining muon tracks must have at least 
Pt > 0 GeV; 

3. Demand exactly one juJ, as defined above. The 
tracks in the fiJ must each have a transverse impact 
parameter |do| between 1 mm and 1.2 m, and the 
distance between the decay vertex and the primary 
vertex must be also between 1 mm and 1.2 m {i.e., 
displaced decay within the tracker); 

4. AR between the prompt muon and each of the 
other muon tracks is less than tt — 0.2 (to veto 
cosmic muons reconstructed as back-to-back muons 
decaying in the tracker); 

5. Veto significant hadronic activity. 

The signal region described above could still be con¬ 
taminated by cosmic-ray (CR) initiated muon bundles 
EH coinciding with a prompt lepton from a pp collision. 
We attempt an estimate of the CR background as follows. 
We use the ATLAS analysis. Ref. EZl, which defines a sig¬ 
nal region as consisting of at least two fiJs in the event. 
In particular, we focus on their validation region, which 
uses a trigger on empty bunch crossings; the events seen 
in that validation sample can only arise from CRs. As¬ 
suming that the probability to find a lepton jet from a 
cosmic ray fact orizes^, the probability to find one CR 
is i^(l/it/) = N 2 observed/-^3/x trigger ~ 6 X 10 and 

the LJ rate is Rlj = P(l/iJ) x Wriggered/Pci?- Here 
Ten ~ (6/15) X (2/3) yr is the approximate livetime 
of the CR data sample acquisition. Since our analysis 
uses a single muon trigger, we rescale Rlj by the ra¬ 
tio of P(l/i trigger)/P(3/i trigger) ^ 180 [78], which gives 
Plj ~ 0.018 Hz. We can now use this probability to es¬ 
timate the number of events in a given livetime period 
T passing a single prompt muon trigger (with pt > 24 
GeV) with a rate R; this prompt muon must coincide 
in a time window AT with a CR-initiated /jJ. We then 
find: 

V(l/i + IjjJ) = RxTx {Rlj x AT) (7) 
« 4 events (8) 

for a 3-year livetime, AT = 25 ns, and assuming an event- 
filter level trigger rate of R = 100 Hz [79] . 

The estimate above should be taken as a heuristic ar¬ 
gument for the CR-initiated backgrounds being a small 
contamination to the signal region. The CR background 


® Due to the fact that a single CR muon is often reconstructed as 
two back-to-back muons at a displaced vertex, it is likely that the 
correlation between finding 1 LJ and 2 LJ is very high, making 
our analysis overly conservative. 
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FIG. 6: 95% confidence level reach of our proposed lepton 
jet and trilepton searches in terms of the sterile neutrino 
simplihed model parameters. The blue lines show the reach 
of the displaced lepton jet search at (dot-dashed) ^/s — 8 
TeV with 20 fb“^, (dashed) ^/s = 13 TeV with 300 fb“^. 
The brown lines show the prompt trilepton reach with (dot- 
dashed) ^/s — 8 TeV with 20 fb“^ and 50% systematic un¬ 
certainty, (dashed) ^/s = 13 TeV with 300 fb“^ and 20% sys¬ 
tematic uncertainty. The thin red dotted line shows the reach 
for the proposed SHiP experiment from Ref. [25]. The shaded 
region is excluded. 


could be further reduced by using the high granularity 
of the tracker and requiring that the two muon tracks 
within the jiJ reconstruct to the same vertex (which was 
not required in Ref. m- Kinematic features such as the 
invariant mass of the jiJ and the alignment of with 
the /iJ could be used to further suppress backgrounds. 
Therefore, we assume a background-free search with in¬ 
tegrated luminosity of 300 fb“^, and define our 2 cr ex¬ 
clusion reach contours by requiring 3 signal events after 
cuts. 

We perform the simulation for the low-mass N signal 
region using Madgraph 5 US]. Because of the all-muon fi¬ 
nal state, we consider only parton-level events. We show 
our estimated sensitivity at the LHC for this signal re¬ 
gion in Fig. both for 8 TeV with 20 fb“^ and for 13 
TeV with 300 fb“^. For masses near ~ 15 GeV, 
the sensitivity of this analysis could be further improved 
by increasing the Ai?o in the definition of {iJ as the N 
decay products’ separation increases. Furthermore, the 
requirement that the fiJ appear at a displaced vertex in 
the tracker {\do\ < Im) could also be relaxed to consider 
DVs in the calorimeters and the muon spectrometer, but 
the background estimate from Ref. has to be modi¬ 
fied for this case. 


IV. PROMPT TRILEPTON SEARCHES FOR 
RH NEUTRINOS 

For masses > 15 GeV, N typically decays 

promptly, and the reconstruction of the decay vertex no 


longer provides significant discriminating power from SM 
backgrounds. In this section, we investigate the most 
promising final states for discovering N in the prompt 
regime. In particular, we find that targeted searches 
in the trilepton final state with no opposite-sign, same- 
flavor (OSSF) leptons can suppress SM backgrounds and 
give a smoking gun signature for lepton-number-violating 
RH neutrinos with Mat ^ While trilepton final 

states have been considered previously for ^ Mw 
and/or Dirac neutrinos (58] |59l | 6 TJ [63U66] . we show 
that the Mat ^ Mw regime presents the LHC experi¬ 
ments with different kinematics than previously consid¬ 
ered. By tailoring the signal selection to the softer kine¬ 
matic regime, we show that trilepton searches have the 
capability of probing Majorana N down to Mat ^ 10 
GeV. 

The only current analysis at the LHC for N in the 
Mat ^ Mw mass range is a CMS search in the 

final state [38|. This search was originally de¬ 
signed for Mat ^ Mw [5QH52] . and has recently been 
re-optimized for Mat ^ Mw [38]. The re-analysis re¬ 
quires two same-sign muons with px > 15 GeV and 
two jets with px > 20 GeV, and seeks to reconstruct 
M^±^±jj Mw‘ It is immediately obvious that, for N 
produced in the decay of , there is insufficient phase 
space to pass all of the kinematic cuts unless the is 
highly boosted; however, if the is boosted, the jets 
in the decay of N are not separately resolved. Therefore, 
this search suffers from extremely tiny signal efficiencies 
for Mat < Mw 0.6 — 0 . 8 %), and for signal events pass¬ 
ing all cuts, one of the jets is typically not from the N 
decay. This can be deduced from the fact that M^±^±jj 
peaks well above Mw for the signal in Ref. [38], whereas 
the correctly reconstructed decay products of N should 
always give a mass below Mw- This suggests that, even 
for signal events, one of the final-state jets is uncorrelated 
with the N decay products, and so the (small) back¬ 
ground looks nearly identical to the signal. Thus, the 
constraints from the search are only comparable 

to or worse than the LEP constraints for Mat ^ Mw- 

Given the challenges in reconstructing both quarks 
from N li^qq' decay as separate jets, we consider in¬ 
stead the purely leptonic decay, M ^ u. 

We propose exploiting the Majorana nature of the ster¬ 
ile neutrino to look for /a^N ^ ja^ji^e^Ue final 

states (see Fig. [^: because there are no OSSF lepton 
pairs in the final state, SM backgrounds involving 7 */^ 
are greatly suppressed. 

Current experimental searches in trilepton final states 
have targeted supersymmetric final states with large 
^x, although CMS has an analysis with low ^x and 
low iLx m- This search has been recast for Mat > Mw 
[65] . and here we recast the analysis to determine the 
constraints on the low-mass signal region Mat ^ Mw- In 
particular, we use the OSSF -0 signal region to find the 
most powerful bound. Using the data from the ^x < 50 
GeV, iLx < 200 GeV, OSSF -0 bin with 0 5-jets from 
Ref. [80], we find that the CMS trilepton analysis is 
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FIG. 7: Production and decay of N for our proposed prompt 
trilepton search with no OSSF lepton pairs. 


competitive with, but does not quite surpass, the LEP 
and CMS analyses for 20 GeV < ^ Mvp. 

Given, however, that a non-optimized analysis in the 
OSSF-0 trilepton channel already gives a bound com¬ 
petitive with other search channels, this suggests that a 
targeted search for N in the trilepton final state would 
give a significant improvement in sensitivity. 

Monte Carlo simulations: We perform a Monte Carlo 
(MC) analysis to estimate the improvement in sensitivity 
that can be obtained with a targeted trilepton search 
for sterile neutrinos. We simulate parton-level processes 
in Madgraph 5 m and shower the events in Pythia 6 
[8T] . For background processes, matrix elements with 
up to two extra partons are simulated and matched to 
the shower using the MLM-based shower-scheme [82] . 
The dominant backgrounds are 7 */Z+ jets, tf, and WZ-\- 
jets. Jets are clustered using the antb/cx algorithm [83] 
implemented with the Fast jet 3 package [84]. Signal 
and background cross sections are normalized to their 
next-to-leading-order values [85H87] . 

A major obstruction to background simulation is that 
the dominant backgrounds for OSSF-0, low-^x and 
\ow-Ht trilepton searches come from processes where 
one or more “non-prompt” (fake) leptons are present 
in the final state. For example, Z/ 7 *+ jets and tt 
backgrounds can fake trilepton signatures if one of the 
final-state jets is mis-tagged as a lepton; this “fake” can 
either come from an actual lepton originating from a 
heavy-fiavor meson decay, or from light hadrons that 
are mis-reconstructed as leptons. Because fake leptons 
are very rare and may rely on improperly modelled jet 
fragmentations, MC estimates for fake leptons are unre¬ 
liable, and the ATLAS and CMS collaborations instead 
use a data-driven approach to estimate lepton fakes 
in their multilepton analyses [38l |80l [88]. Since we do 
not have access to the resources needed for data-driven 
estimates, we adopt an approach proposed by Ref. [89] . 
which takes jet-enriched samples and constructs a map 
from jet kinematics to fake lepton kinematics. This 
method allows for the use of reasonably sized samples of 
Z/ 7 *+ jets and tt events to obtain sufficient statistics 
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GeV, < 50 GeV bin. For both signal mass points, there 
is a cutoff in the distribution at Mw, with the peak more 
prominent for smaller Mn- 


for estimating fake lepton backgrounds. We describe the 
procedure and validation of this method in Appendix [Aj 

Signal kinematics: We apply basic selection criteria 
similar to the OSSF-0, 0 6-jet bin for the CMS trilepton 
analysis [80]. Requiring a leading lepton with px > 20 
GeV and all subleading leptons with px > 10 GeV, we 
demand exactly three leptons, zero OSSF lepton pairs, 
and zero 6-tagged jets with px > 30 GeV (using the b- 
tagging working point from Ref. [80]). Defining iLx as the 
scalar px sum of jets with px > 20 GeV, we apply upper 
cuts on iLx and ^x to suppress tt backgrounds. For the 
histograms shown below, we apply Ht < 200 GeV and 
^x < 50 GeV, although stricter cuts are applied for the 
final analysis. 

Taking as our convention that if (if) is the hardest 
(softest) same-sign lepton, and i^ is the lepton of op¬ 
posite sign, we study various kinematic distributions of 
the charged leptons. In particular, we find two observ¬ 
ables that are powerful discriminants between signal and 
background. The first is the trilepton invariant mass, 

; because the invariant mass of the three leptons 
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FIG. 9: Histograms of in the OSSF-0, 0-6, Ht < 200 

GeV,^T < 50 GeV, < 80 GeV bin. For both signal 

mass points, the signal peaks below Mat, with the peak more 
prominent for smaller Mat. 


plus the neutrino reconstructs the TV, this distribution 
has a sharp cutoff at Mvp, as we show in Fig.|^. Indeed, 
for masses Mat ^ Mw^ most of the energy from the W 
decay must go into the charged leptons for a high effi¬ 
ciency of passing kinematic cuts, and therefore a strong 
peak is observed even though one of the final states is in¬ 
visible. For Mat ^ Mw, all of the leptons (including the 
neutrino) arising from the N decay are relatively hard, 
while the muon from the original W decay is soft in turn, 
and so the W peak is not as pronounced. 

The second useful discriminant is the invariant mass 
To clearly show its effects, we impose a cut of 
< 80 GeV in addition to the selection criteria 
described above. For Mat ^ Tfw, if and i^ both come 


^ In the CMS search for semileptonic decays of V —)■ fj,^ [38], 
this type of kinematic observable is used, but because one of the 
jets is typically from uncorrelated radiation in signal events, the 
signal peaks above M\y and the observable is not an effective 
discriminant of signal from background. This is in contrast with 
our findings for trilepton searches. 


from the decay of N, and so their invariant mass is 
kinematically restricted to be < Mat, which we show 
in Fig. For Mat ^ Tfw, a peak structure is again 
obtained, while for Mat ^ all three charged leptons 
have comparable momenta, and it is no longer possible 
to identify which lepton came from the N decay. Thus, 
the distribution is broader and looks more like 

background, degrading the discriminating power (as 
shown in Fig. [^with Mat = 50 GeV). 

Proposed analysis: Taking into account all of the 
above, we impose the following kinematic selections: 

1. Two prompt isolated same-sign muons and one 
prompt isolated opposite-sign electron with px > 
10 GeV (leading px > 20 GeV); 

2. Ht < 50 GeV and ^x < 40 GeV (to further sup¬ 
press fake tt and Z tt backgrounds)^; 

3. < 80 GeV (for Mat < 50 GeV, we addi¬ 
tionally require > 60 GeV); 

4. An upper cut on separately optimized for 

each value of Mat. 

We optimize the cuts on for each Mat to maximize 

signal significance for ^/s = 13 TeV and 300 fb“^, scan¬ 
ning over values of the cut ranging from 10-60 GeV. The 
optimal cut is typically just below Mat- We determine 
the signal reach at 95% confidence level, with the num¬ 
ber of signal events exceeding max(3, 2crbkd) foi* a given 
integrated luminosity. We include both statistical uncer¬ 
tainty and a projected systematic uncertainty in abkd- 

The projected reach is determined both for = 8 
TeV and 20 fb“^ (50% systematic), as well as ^/s = 13 
TeV and 300 fb“^ (20% systematic)^. For comparison, 
the current CMS search for Majorana neutrinos in the 
2/i^ + 2j channel is dominated by fakes and the system¬ 
atic uncertainty on the background estimate is ^ 30% 
j38] . We show our estimate of the trilepton search reach 
in Fig. For 10 GeV < Mat < Mw, the prompt 
trilepton search can substantially improve upon the LEP 
bound, and a search with Run I data may already cover 
new territory. The sensitivity diminishes at low mass 
because the leptons from N decays originate from a dis¬ 
placed vertex; this region is best covered by the lepton 


® Our simulations indicate that, for signal, and^.^ may peak 
even lower than this; however, resolution effects, initial state ra¬ 
diation, and pile-up become important for smaller values of these 
variables, and so the optimal cut value should be determined by 
ATLAS and CMS. 

^ In a data-driven background estimate, an increase in statistics 
in the control region will lead to a smaller uncertainty in the 
extrapolation into the signal region, which explains our choice 
of a smaller systematic uncertainty with the larger dataset at 13 
TeV. 
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jet search proposed in Section |III| In the overlapping 
region, performing both searches can improve sensitiv¬ 
ity to N and increase confidence in a putative signal. 
For Mat ^ Mw: the reach weakens because the decay 
jj^^N becomes phase-space-suppressed. 

In summary, the search for N in prompt trilepton fi¬ 
nal states nicely complements the proposals for long-lived 
searches at fixed-target experiments and at the LHC, and 
is a promising channel for probing RH neutrinos in both 
Run I and Run II data. 


V. DISCUSSION AND CONCLUSIONS 

Theories with right-handed neutrinos at or below the 
weak scale can resolve many of the shortcomings of the 
SM, providing an explanation for neutrino masses, the 
origin of the baryon asymmetry, and a candidate for dark 
matter. We have shown that, in the well-motivated mass 
range rric < Mn < Mw , the LHC can serve as a powerful 
probe of right-handed neutrinos due to the large number 
of W bosons produced, complementing the lower-mass 
searches at fixed-target probes such as the proposed SHIP 
experiment. In particular, the leptonic interactions of 
a Majorana neutrino give rise to distinctive multilepton 
searches in displaced and/or prompt final states. 

We have proposed two search strategies that can im¬ 
prove the reach of right-handed neutrinos by up to two 
orders of magnitude in cross section: a search for a dis¬ 
placed lepton jet in association with a prompt lepton, 
which covers the mass range rric < ^ 20 GeV, and 

a prompt trilepton search that covers the mass range 
10 GeV < Mat ^ Mw Optimizing the kinematic se¬ 
lections with respect to background can lead to clean 
signatures with good discovery potential, even though 
the decay products of right-handed neutrinos are rel¬ 
atively soft. Because the final states are entirely lep¬ 
tonic, the dependence on pile-up and other factors with 
high-energy/luminosity running should not significantly 
degrade the prospects for the discovery of right-handed 
neutrinos. 

Before concluding, we discuss briefly the prospects for 
right-handed neutrino searches in the range Mat ^ 
at ATLAS/CMS. Here, the dominant production mech¬ 
anism is not from W decays, but from H-hadron decays 
such as D/i+V. Naively, it seems that the strong 

production cross section of N in this mass range could 
lead to reasonable discovery prospects. Furthermore, the 
RH neutrinos are typically very long-lived in this mass 
range (cr > 1 — 10 m), and so one can take advantage of 
the large ATLAS/CMS detectors to look for N decays in 
the muon spectrometer m- However, we find that for 
^6 — '^c ^ Mn ^ ^ 6 , the tat results in decays inside the 
muon spectrometer, but the B is kinematically forced to 
decay through its mixing with the up quark, and this 
leads to a suppressed N production rate. Conversely, 
for Mat ^ ruiy — rUc, the branching fraction of H ^ V is 
larger, but the lifetime is longer so that one must pay a 


severe penalty in rate to require the decay length to be 
< 10 m. Combining this with the requirement of having 
sufficiently hard final states to reconstruct the lepton 
from B decay and the displaced vertex inside the muon 
spectrometer, we estimate the projected sensitivity to be 
worse than the current constraints. Thus, the lepton-jet 
final state from W decay is still the best way to search 
for RH neutrinos for Mat ^ mij with ATLAS/CMS, 
although direct searches of N production from B decays 
at H-factories and LHCb will have sensitivity in this 
mass range [36l |9T1 [92] . 
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Appendix A: Fake Lepton Simulation 


We use the method of fake lepton simulation proposed 
in Ref. [89| . This approach relies on the fact that fake lep¬ 
tons are initiated by jets, and the leptons consequently 
inherit the kinematic properties of the jet. Therefore, 
fake lepton kinematics can be estimated from jet-rich MC 
samples by mapping jet kinematics onto the fake lepton 
momenta. In particular, the authors of Ref. [89] found 
that lepton fakes can be well-estimated by applying: a 
mistag rate^ which is a function of the jet px; 

and a transfer function^ which is a probability dis¬ 

tribution function that maps pxj into the fake px^- These 
have the functional forms: 


^j^iiPTj) — ^200 


Tj^eia) = ^exp 


(1 - rio) 
20-2 


200-pTi/GeV 


200 - 10 


(M) 

(A 2 ) 


In this expression, px^ = (1 — (a)pxj, ^200 parameterizes 
the overall mistag rate, rio gives the gradient for the px- 
dependent mistag rate, and p and a are the mean and 
standard deviation of the Gaussian distribution mapping 
the jet momentum to the lepton momentum. If a 7 ^ 0, 
some of the jet momentum is not mapped into the lepton 
momentum, and it is a ssum ed to contribute to ^x- The 
Gaussian function Eq. (A 2 ) is truncated so that 0 < a < 
1 , and AT is a normalization factor. 

We find that agreement with experimental data is ob¬ 
tained only if we assume that all fake leptons arise from 
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FIG. 10: Comparison of our MC simulation of 7 */Z+ jets 
(top, red), tt (middle, green), and IFZ+jets (bottom, blue) 
with the observed values from the CMS trilepton analysis 
(black dots) [80]. The analysis is at ^/s — 8 TeV, £ = 20 fb“^, 
and applies all of the selections in Ref. m including the cuts 
< 50 GeV and Ht < 200 GeV. The bins are divided ac¬ 
cording to the number of 5-tagged jets, whether there exists 
an OSSF lepton pair, and the mass of any OSSF lepton pair. 
The bins are: 1) 0-5, OSSF-1, M^+^- < 75 GeV; 2) 0-5, OSSF- 
1, \M^+^- - Mz\ < 15 GeV; 3) 0-5, OSSF-1, M^+^- > 105 
GeV; 4) 0-5, OSSF-0. Bins 5-8 are the same as bins 1-4, but 
with at least one 5-tagged jet. 



1 2 3 4 5 6 7 8 

bin number 


jets containing heavy-flavor hadrons. Consistently, we 
find that a good fit to data is obtained with the inputs 
rio = 1, /i = 0.5, and a = 0.3. These parameters give a 
flat mistag rate in jet px, and equally divide the “jet” en¬ 
ergy between the fake lepton and the neutrino, which is 
expected from leptons originating in heavy-flavor decays. 
These parameters are somewhat different from those ob¬ 
tained in Ref. [89| , but in both their analysis and our own, 
we find that changing the parameters in the efficiency and 
transfer function does not substantially change the fit, 
provided the total normalization remains fixed. This is 
due to approximate degeneracies present among the var¬ 
ious parameterizations of the fake lepton functions: for 
example, a softer momentum spectrum can be obtained 


both with a large /i, in which case the lepton momentum 
is imparted with less of the original jet momentum, or 
rio ~ 1 , which gives a larger weighting to soft momenta 
relative to rio = 0. In our comparison with the data 
described below, we find that varying the fake simulator 
input parameters changes the contributions to various 
bins by a factor < 2 , and similarly, the reach in |V^ArP 
of our trilepton analysis similarly changes by < 2 . 


As our paper focuses on trilepton signatures, we fix 
the overall fake efficiency using the CMS trilepton data 
[80] : we normalize our background estimates to the 
OSSF- 1 , 0 5-jet, M£+£- ^ Mz bin (“bin 2 ” in Fig. 10), 
and find 6200 = 4.6 x 10“^. Having fixed all of the pa¬ 
rameters in the fake lepton simulator using data from a 
single bin, we can extrapolate the results to all bins, and 
we show the result in Fig. We see that the combi¬ 
nation of prompt and fake lepton simulated backgrounds 
matches extremely well the backgrounds from Ref. m- 
If we had not made the assumption that fake leptons 
originated exclusively from heavy-flavor jets, the relative 
rates of 7 */Z+ jets and tt would change by an order of 
magnitude, and the agreement would no longer be ac¬ 
ceptable. 


Although not shown here, we also validated our fake 
lepton simulation against other ATLAS and CMS same- 
sign lepton searches |38l|88]. In these analyses, only two 
final-state leptons were required (both of the same sign), 
and consequently contributions from multijet processes 
with two fake leptons are also important. The fake sim¬ 
ulation procedure described above also works well at re¬ 
producing the spectra from Refs. [38l|88] (except, on oc¬ 
casion, at the very high-px part of the background spec¬ 
trum), but the overall normalization typically must be 
adjusted for each analysis. This is not surprising, since 
each analysis looks at a different final state (Ref. [ 88 | is 
an inclusive same-sign lepton search, while Ref. [38] re¬ 
quires two additional jets). It does suggest that the good 
shape agreement between MC simulation and CMS data 
shown in Fig. is robust provided we do not change the 
basic final-state topology considered in Ref. [80], and so 
we preserve the basic object selection from that analysis. 
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